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I.  Introduction 


•''-In  distributed  tactical  decision  making,  commanders  make  decisions  based  on  data  gen¬ 
erated  by  their  local  territory  as  well  as  on  decisions  and  information  communicated  to  them  by 
other  commanders.  The  timely  aspect  of  the  decisions  is  controlled  by  the  rate  with  which  the 
commanders  receive  data,  which  is  in  turn  controlled  by  the  deployed  transmission  algorithms, 
(such  as  those  in  [3]  and  (4)).  The  accuracy  of  the  decisions,  on  the  other  hand,  depends  heavily 
on  the  way  that  the  decision  makers  perceive  their  environment.  This  perception  corresponds  to 
a  number  of  alternative  models,  where  the  latter  are  a  priori  developed  based  on  the  various  bat¬ 
tle  scenaria.  Each  such  model  is  associated  with  an  appropriate  decision  mode.  Furthermore,  a 
model  (or  battle  scenario)  may  shift  to  another  such  model,  at  a  random  time.  In  the  event  of 
such  a  shift,  it  is  thus  crucial  that  the  commander  be  alerted,  for  adaptation  to  the  appropriate 
decision  mode.  In  this  paper,  we  describe  an  algorithm  which  can  be  deployed  in  an  automatic 
fashion  ,  to  alert  the  commanders  about  changes  in  battle  scenaria. 


II.  The  Algorithm 

Each  battle  scenario  is  modelled  by  a  stochastic  process.  Let  the  various  battle  scenaria  be 
indexed,  say  from  1  to  M.  Then,  the  ith  battle  scenario  is  represented  by  a  distribution 
Fjfx!,..,  xn,..),  where  x,,..,  xn,...  are  the  data  outputs  observed.  Let  us  assume  that  the  above 
distributions  possess  density  functions,  and  let  then  fi(x„  lxi,..,xn_,)  denote  the  density  of  the 
datum  xn,  given  the  data  xj ...  xn_, ,  and  given  that  the  ith  battle  scenario  is  present. 

Let  the  jth  battle  scenario  be  known  to  be  present  at  time  zero.  Our  objective  is  to  decide  if 
a  change  in  this  scenario  occurs,  and  if  the  answer  is  yes,  then  to  decide  what  is  the  new 
scenario.  Furthermore,  if  a  change  does  occur,  we  wish  to  decide  affirmatively  as  soon  from  the 
instant  of  its  occurence  as  possible.  Let  us  define. 


g.j  (xi...,  xn)  =  log 


fj  (xn  I  X|  , ..,  Xn_| ) 
fj(xn  I  X, ,..,  Xn_j  ) 


Tfj1'  such  that. 


T{jn)  =  0 


T|"+1)  =  max(0,  T}]”  +  g^x,,....  xnH), 

A  set  (5,j  ,  l<i<\l)  of  positive  thresholds,  where  5„  is  associated  with  the  sequence 
{T<n).n>l} 

The  thresholds  {5()  ,l<i<M)  are  determined  a  priori  based  on  false  alarm  and  correct  deci¬ 
sion  criteria.  The  algorithmic  decision  is  then  as  follows: 

Decide  that  scenario  j  shifted  to  scenario  k,  the  first  time  n,  such  that  T^'  >  5^,  while  for 
each  i  different  than  k,  the  value  Tjj1>  is  still  below  its  threshold  8ir 
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The  above  algorithm  is  dearly  sequential,  and  numerically  very  efficient.  Furthermore,  as 
proven  in  fl],  it  is  also  asymptotically  optimal;  that  is,  there  is  no  other  algorithm  that  can  detect 
a  change  in  model  faster  than  the  above  algorithm.  In  addition,  if  there  is  lack  of  confidence  in 
some  of  the  observed  data  or  ambiguity  in  the  models  of  the  battle  sccnaria,  then  the  algorithm 
can  be  easitly  modified,  to  basically  maintain  its  performance  characteristics,  (see  [2]). 


III.  Example 


Let  the  observed  data  be  binary  and  independent,  and  let  us  consider  two  scenaria,  I  and  II. 
Let  :L  be  known  that  under  scenario  I,  a  bit  has  value  1  with  probability  p,  and  has  value  0  with 
pro  ,uodity  1-p.  Let  it  also  be  known  that  under  scenario  II,  a  bit  has  value  1  with  probability  q, 
and  has  value  0  with  probability  1-q.  Let  scenario  I  be  initially  present,  and  let  X),  x2, . . . , 
denote  the  observed  binary  sequence.  Then,  for  detecting  a  change  to  scenario  II,  the  algorithm 
easily  reduces  as  follows: 

Select  some  threshold  5  positive. 


Set  T0  =  0  and  Tn  =  max(0,Tn_i  +  xn  + 


log 


i-q 

1~P 


log 


qn-p) 

p(l-q) 


) 


Decide  that  scenario  I  shifted  to  scenario  II,  the  first  time  n  such  that  Tn  >  6. 

Here,  we  will  study  the  nonasymptotic  performance  of  the  algorithm.  In  particular,  given 
the  threshold,  S,  we  will  anlayze  the  false  alarm  and  the  power  characteristics  induced  by  the 
algorithm.  Let  us  define. 


Pr,5(n):  The  probability  that  the  threshold,  5,  is  first  crossed,  at  the  time  instant,  n, 

under  the  condition  that  the  p  to  q  change  actually  occurred  ( 4 ) 

just  after  the  datum  xr. 


The  probabilities,  |P-,6(n) ;  n  >  1 },  represent  then  a  false  alarm  set,  while  the  probabilities, 
(Po.5(n) ;  n  >  1 },  represent  a  power  set.  Indeed,  the  former  probabilities  show  how  probable  is 
the  crossing  of  the  threshold,  at  different  time  instants,  n,  when  the  p  to  q  change  never  occurs. 
The  latter  probabilities  correspond,  instead,  to  threshold  crossings,  at  dilfermt  time  instants,  n. 
when  the  "alarming",  parameter-q,  process  is  active,  since  after  time  zero.  We  note  that  when 
the  threshold,  5,  is  finite,  it  is  eventually  crossed  with  probability  one,  even  when  the  p  to  q 
change  never  occurs.  The  false  alarm  curves,  {P-,5  (n);n>l ),  and  the  power  curves. 
(Po.5  (n) ;  n  >  1 ),  can  be  only  computed  numerically.  For  efficiency  in  such  computations,  we 

i  1  -q 
°B  Tm 

apply  an  approximation.  In  particular,  we  approximate  the  constant, - - —  by  the  ratio  of 


two  integers.  That  is,  we  write. 


(5) 


C 


C 


c 


log 


JlSL 

1-P 


log 


qd-p) 

p(i-q) 


l 

s 


;  where  /  <  s,  and  /  and  s  are  both  natural  numbers. 

Then,  we  write  the  algorithm  in  the  following  equivalent  form,  where,  without  lack  in  gen¬ 
erality,  the  threshold,  t,  in  the  latter  is  a  natural  number. 

Select  the  threshold,  t.  Decide  the  change  has  occurred,  at  the  first  time,  n,  such  that: 
Tn  >  t,  where: 

T0  =  0,  Tn=max(0,  Tn_) +y„)  (6) 

yn  =  s  xn  -  /  (7) 


Clearly,  the  updating  step,  yn,  in  (7)  is  an  integer.  It  is  equal  to  -/,  when  xn  =  0,  and  it  is  equal 
to  s  -  /  >  0,  when  xn  =  1.  Let  us  denote  by,  Pril(n),  the  probability  in  (4),  as  applied  to  the  modi¬ 
fied  algorithm  in  (6)  and  (7).  It  is  then  possible  to  obtain  recursive  expressions  for  the  computa¬ 
tion  of  the  probabilities,  {Pril(n) ;  n  £  1 },  via  a  Markov  chain  formulation.  The  key  element  in 
the  latter  formulation  is  the  probability,  Pr,t(n,j),  that  T'(xk)  <t;¥k<n  and  T(xn)  =  j,  given 
that  the  p  to  q  change  occurs  just  after  the  datum,  xr.  Indeed,  we  find. 


i.  If,  t  -  1  >s>/+  I: 

Pr,t(n.O)  =  ( 1  — v)  2  PM(n-l,  i) 

i=0 


Pr,t(n.j)  =  (1-v)  Pr.,(n-1,  j+/) ;  1  <  j  <s  -/ 


(8) 


Pr.i(n.j)  =  (1-v)  Pr.c(n-1,  j+/)  +  V  Pr,t(n-1,  j-s  +  /);  s-/  <j  £  t— 1— / 
P:.i(n,j)  =  v  PM(n-l,  j-s+o  ;  t— 1— /  <  j  <  t  -  1 


ii.  If,  t  -  1  >  s  =  /+  1: 


Pr .i(n,0)  =  (1-v)  £  Pr,i(n-1,  i) 

i-0  (9) 


Pr.t(n.j)  =  (l-v)Pr,l(n-l.  j+/)  +  vPr  l(n-l,  j-1);  1  <  j  <  t  -  I  —  / 


s 

Pr.i(n.j)  =  V  PMU1-1,  j-1)  it  —  1  — /<j<t  —  1  ! 


c 


I 
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;  where. 


v 


p  ;  if  n  £  r 
q  ;  if  n  >  r 


(10) 


Pr.i(n)  =  v  £  Pr.i(n-l,j);t-l>s>/+l  (11) 

j=i  -s+l 

The  recursive  expressions  in  (8)  and  (9)  allow  for  the  efficient  computation  of  the  probabil¬ 
ity,  PM(n)  in  (1 1),  for  every  r,  and  every  n. 

The  false  alarm  probabilities,  {P«,t(n) ;  n  >  1 },  and  the  power  probabilities, 
{P0>t(n) ;  n>  1),  evolve  as  a  special  case  of  the  probabilities,  {Pr,t(n)),  and  they  can  be  also 
computed  recursively  then.  Let  us  now  define  the  following  two  quantities. 

r. 

a(p,q,n,t)  =  £  P„.,(i)  (12) 

i=0 

n 

P(p,q,n,t)  =  £  P0.t(i)  (13) 

i=0 

The  quantities,  a(p,q,n,t)  and  p(p,q,n,t),  are  the  probabilities  that  the  threshold,  t,  is 
crossed  by  the  algorithm  in  (6)  and  (7),  before  or  at  the  time  instant,  n,  given  that  all  the  data  are 
respectively  generated  by  the  p-parameter  and  the  q-parameter  process.  Thus,  a(p,q,n,t)  and 
P(p,q,n,t)  are  respectively  the  false  alarm  curve  and  the  power  curve  induced  by  the  algorithm  in 
(6)  and  (7).  Those  two  curves  represent  the  nonasymptotic  performance  of  the  algorithm.  We 
computed  those  curves  for  p  =  .01  and  q  =  .05,  and  for  various  choices  of  the  threshold,  8.  In  our 
computations,  we  applied  the  rational  approximations  represented  by  expression  (5)  and  the 
algorithm  in  (6)  and  (7),  and  we  used  the  recursions  in  (8)  ad  (9).  Our  results  are  shown  in  Fig¬ 
ure  1.  From  this  figure,  the  powerful  nonasymptotic  characteristics  of  the  algorithm  are  evident. 
For  relatively  small  sample  sizes,  it  simultaneously  attains  both  low  false  alarm  and  high  power. 
The  choice  of  the  threshold,  8,  will  be  determined  from  the  specific  false  alarm  and  power 
requirements,  and  it  will  depend  on  where  the  emphasis  (power  versus  false  alarm)  is  placed. 
When  the  probabilities,  p  and  q,  arc  further  apart  from  each  other,  than  the  pair,  p  =  0.1  and  q  = 
.05,  then  the  false  alarm  and  the  power  curves  are  drawn  further  apart,  than  those  in  Figure  1. 
The  performance  characteristics  of  the  algorithm  are  then  further  improved.  Studies  similar  to 
those  in  this  example  can  be  performed,  when  dependent  binary  data  structures  are  instead  con¬ 
sidered,  (see  (51),  or  when  Gaussian  data  are  present. 


IV.  CONCLUSIONS 

In  this  paper,  we  exhibited  a  simple  sequential  algorithm  which  we  have  previously  fully 
analyzed,  (see  [  1  j  and  [2]),  and  which  can  be  deployed  to  detect  changes  in  battle  scenaria.  The 
algorithm  can  be  implemented  in  an  automated  fusion.  At  the  same  time,  it  would  be  interesting 
to  see  how  closely  can  trained  human  subjects  approximate  the  algorithmic  operations  and  how 
fast  can  they  react  to  the  algorithmic  alerts.  We  propose  that  experiments  along  those  lines  be 
performed  by  the  researchers  on  human  factors,  within  the  DTDM  program. 
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degrees.  In  addition,  courses  in  the  humanities  are  offered  within  the  School. 

The  University  of  Virginia  (which  includes  approximately  2,000  faculty  and  a  total  of  full-time 
student  enrollment  of  about  1 6,400),  also  offers  professional  degrees  under  the  schools  of  Architecture, 
Law,  Medicine,  Nursing,  Commerce,  Business  Administration,  and  Education.  In  addition,  the  College 
of  Arts  and  Sciences  houses  departments  of  Mathematics,  Physics,  Chemistry  and  others  relevant 
to  the  engineering  research  program.  The  School  of  Engineering  and  Applied  Science  is  an  integral 
part  of  this  University  community  which  provides  opportunities  for  interdisciplinary  work  in  pursuit 
of  the  basic  goals  of  education,  research,  and  public  service. 
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